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ABSTRACT 

We present the first spatially resolved spectroscopic imaging observations 
of the 12 C0 (1-0) line emission in the central 2.5 kpc of the Seyfert 1 galaxy 
NGC 4151, obtained with the IRAM Plateau de Bure Interferometer (PdBI). 
Most of the cold molecular gas is distributed along two curved gas lanes about 
1 kpc north and south of the active nucleus, coincident with the circumnuclear 
dust ring noted by previous authors. These CO arcs lie within the Inner Lindblad 
Resonance of the large scale oval bar and have kinematics consistent with those 
derived from neutral hydrogen observations of the disk and bar. Two additional 
gas clumps are detected that show non-circular motions - one associated with 
the southern gas lane and one lying ~600 pc north of the nucleus. Closer to the 
nucleus, no cold molecular gas is detected in the central 300 pc where abundant 
near-IR H2 line emission arises. This suggests that the H2 line emission is not 
a good indicator for a cold gas reservoir in NGC 4151 and that the H 2 is likely 
photo-excited by the AGN. The upper limit of the CO mass in the central 300 pc 
is sufficient to support the AGN activity at its current level for 10 7 yrs. The total 
cold molecular mass detected by PdBI is 4.3 x 1O 7 M . Finally, 3 mm contin- 
uum emission arising from the location of the AGN is detected with a flux of 
Szmm ~ 14 mJy and appears to be unresolved at an angular resolution of 2'.'8 
(~180 pc). 

Subject headings: galaxies: Seyfert - galaxies: ISM - galaxies: kinematics and 
dynamics - galaxies: individual(NGC4151) 
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Introduction 



Nuclear activity in galaxies is thought to be driven by the release of gravitational poten- 
tial energy from material accreted by a central supermassive black hole. In this scenario, fuel 
transportation mechanisms in the nuclear regions must be effici ent on time scales comparable 



to th at of the nuclear activity of the order of 10 7 -10 8 yrs (e.g. iMartinil 120041 : iMarconi et al. 



20041 ) . Stellar b ars are well known to drive gas from the outer large-scale disk to the cen 



tral kiloparsec ( Sakamoto et al.l Il999l ; ISheth et al.l 120051 ) and the bar-within-bar scenario 



( jShlosman et al.l Il989l ) was proposed as a mechanism to deliver gas closer to the central 
black hole. Recently, other gravitational perturbations, such as m=l modes or gas density 
waves, have been suggested to play a major role in the transport of gas to the most c entral 
regions (e.g. 



Englmaier fc Shlosmanll2000l ; lEmseilemll200ll ; iGarcia-Burillo et al.ll2003af ) 



Molecular gas dominates the interstellar medium (ISM) in the centers of nearby spiral 
galaxies, therefore CO lines are the best tracers for the nuclear gas distribution and dynamics. 
High resolution CO observations of the central kpc regions of nearby active galaxies provide 
insights on the ISM properties in these regions and the processes that could move gas inwards. 
The recently complete d NUGA survey with the IR AM mm- interferometer PdBI (Plateau de 



Bure interferometer) (IGarcia-Burillo et al.l l2003al ) observed 12 Seyfert and LINER galax- 



ies with different nuclear and host properties and showed that molecular gas is present in 
the centers, however with diverse morphology and kinematics. This already suggests that 



ing the blac 


t hole ( 


Garcia-Burillo et al. 


2003a. 


20051 


Com 


)es et al. 


2004 




Krips et al. 


2005; 


Boone et al. 


2007; 


B 


!unt et al.l 


2008; 


Lindt-Kriee et al. 


2008 


Casasola et al. 


2008: 


Haan et al. 


2009; 


Combes et al. 


2009 


). Thus, the mechanism for transporting material from the central 



few hundred pc to the SMBH in order to induce and sustain the nuclear activity remains 
debated. Observations of molecular gas in nearby AGN can begin to bridge this gap, con- 
necting large and small spatial scales and examining the next link in the fueling chain. In this 
context, the nearby spiral galaxy NGC4151 hosting a Seyfert 1 nucleus is an ideal target. 



NGC4151 is one of the best studied Active Galactic Nuclei (AGN) (see r eview by lUlrich 



20001 ) and situated at a distance of 13.3M pc (1" = 65pc. iMunde 



spiral galaxy, classified as (R')SAB(rs)ab ( Ide Vaucouleurs et al. 
on (i=21°). In addition to a remarkably large stellar bar (3'x2' 



1 et al.ll2003l ). It is a barred 
199ll ) , and seen nearly face- 
11.7kpc x 7.8 kpc), two 



faint spiral arms extend out to a radius of ~6' (23.4 kpc). It has been extensively observed at 
all wavelengths and on various spatial scales. In particular neutral hydrogen (H i) traces the 



1 Based on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported 
by INSU/CNRS (France), MPG (Germany) and IGN (Spain). 
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large-scale stellar bar and the outer spiral arms. On kpc scales, IVila-Vilaro et al.l (119951 ) have 
detected two red arc-like features, possibly produced by dust extinction of the background 
stellar continuum, which delineate a circumnuclear elliptical ring (hereafter: central dust 
ellipse) with a semi-major axis of ~18" (~1.2kpc). In addition a large Extended Narrow Line 
Region (ENLR) extends up to 20" (1.3kpc) and is consistent with ambi ent galactic gas being 
photo- ionized by a cone of nuclear UV radiation ( iPenston et al.lll990f). Evidence for radial 



outflow have been observed and modelled in the NLR of NGC4 151 (Heck man fc Balick 



Schulzlll990l;lHutchings et al. 



1998; Kaiser et al. 2000: Crenshaw et al. 



2000 



Das et al. 



1983 



2005 



Storchi-Bergmann et al.ll2010l ). The bright Seyfert nucleus of type 1.5 displays rapid optical 



varia bility, and a small ~ 600 pc long radio jet at a position angle (PA) of ~ 80°( lPedlar et al. 



19931 ). On scales of about 1" to 2", optical CCD imaging revealed evidence for obscuration in 
the form of a reddened band of enhanced extinction (or low ionization) crossin g the central 



regio n , and aligned approx imately perpendicular to the ENLR and radio jet (IPerez et al. 



19891 : iTerlevich et al.l Il99ll ). In addition, there i s also evidence of neutral gas and dust 



within ~10pc of the nucleus ( IMundell et al.lll995l . 119991 1 . A summary of the properties of 
NGC4151 is provided in Tab. U 

Here we present the first mm-interferometric 12 CO(1-0) observations of the central 2 kpc 
of NGC4151, encompassing the dust ellipse, obtained with PdBI. These data provide a link 



betw een the large-scale bar and the (putative) torus seen in H I in absorption (IMundell et al. 



20031 ) . We describe the observations and data reduction in Sect. [2j The resulting CO maps 
and molecular gas properties are presented in Sect. [3J The relation between the molecular 
gas and the large-scale properties of NGC4151 is analyzed in Sect. 14.11 In Sect. 14.2} cold 
and warm molecular gas properties in the innermost region are compared and we discuss our 
results in the context of AGN fueling in Sect. 14.31 Finally we conclude in Sect]5j 



2. Observations and Data Reduction 

The 12 CO(1-0) emission line at 3mm was observed with IRAM/PdBI using 6 antennas 
in C and D configuration on December 26 and 27, 2008 and April 9 and 19, 2009, respec- 
tively, providing baselines from 24 m to 176 m. The calibration and map ping were done with 



the s tandard IRAM GILDAS software packages CLIC and MAPPING ( IGuilloteau fc Lucas 



20001). The phase center of the observation was at RA= 12 h 10 m 32. s 579 and DEC= +39 o 24'20"63 



-i 



in the J2000.0 coordinate system, all velocities are observed relative to v sys = 995 kms 
The quasars 3C273, 3C84, 3C345, 1144+402, 1308+326, 0923+392, 0528+134 and 1156+295 
were observed as flux calibrators. During all observations, 1144+402 was used as phase cal- 
ibrator and observed every 24min. 
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We separately calibrated each day of observations using the corresponding phase and 
flux calibrators. All observations were then merged into a single dataset. The correlator was 
set to a spectral resolution of 1.25MHz per channel, corresponding to a velocity resolution 
of 3.25 km s -1 . We binned by two channels for a final spectral resolution of 6.5 km s™ 1 . In 
addition, a pure continuum map was created by averaging all channels with no line emission. 
In order to obtain a data cube containing only emission line information, the continuum was 
subtracted from the line data in the uv-plane using the task UV_SUBTRACT in GILDAS. 
The continuum and emission line only data were cleaned separately, using uniform and 
natural weighting (see Tab. [2]). All cleaned data have a pixel size of 0.5", with an image size 
of 256x256 pixels . 

The channel maps of the naturally weighted CO (1-0) line are presented in Fig. [1] along 
with the dirty beam and the uv coverage. All data are presented without primary-beam cor- 
rection applied. At the observed frequency, the typical FWHM of the PdBI primary beam is 
44". Hereafter, we present the continuum emission at the uniform resolution (2.83" x2. 16") 
and the results from the naturally weighted emission line cube at 6.5 km s -1 velocity res- 
olution, which has a lower spatial resolution (beam of 3.44" x 2.96" or 224pcxl92pc) but a 
better sensitivity, with a noise level of 2 mJybeam -1 , than the uniformly weighted cube (for 
which the rms is 2.8 mJybeam -1 ). 



3. Continuum, Molecular Gas Distribution and Kinematics 



Here we present the properties of the mm continuum and CO (1-0) line emission in the 
central two kpc of NGC4151, as observed with PdBI. 



3.1. Millimeter continuum emission 



The cleaned millimeter continuum map derived with uniform weighting is shown in 
Fig. [2j The beam size is 2.8"x2.2" and the rms noise is 0.3 mJybeam -1 . The total inte- 
grated flux density is 14 mJy and the peak flux is 12 mJybeam -1 . The millimeter continuum 
peaks at RA= 12 fe 10 m 32 s 56 a nd DEC= +39°24 / 21'. / 06 (J2000.0) corresponding to the AGN 
location ( IMundell et al.ll2003l ). with an offset of 0.25", consistent with our astrometric un- 
certainty of about 0.4". At this resolution it is barely resolved: a two-dimensional Gaussian 
fit gives a convolved source size of 3.2" x 2.3" with a PA of 93°. While the elongation of the 
continuum is similar to the shape of the clean beam, this suggests a small extent in east-west 
direction roughly consistent with the orientation of the radio jet that is elongated along a 
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PA of 77° (IMundell et al.lll995[ ) 



3.2. CO line emission 

We used the task MOMENTS in GIPSY to construct the intensity, velocity and velocity 
dispersion maps of the CO(l-O) line emission, using a flux clipping level of 3er=6 mJybeam -1 
per channel and considering emission only as real if it appears at least in two consecutive 
channels. The width of 2 channels is 13kms _1 , smaller than the typical line width (20kms _1 ) 
seen throughout the data cube. The three moment maps of the naturally weighted data cube 
are shown in Fig. [3J 



3.2.1. Molecular gas morphology and masses 

The integrated intensity map of the CO line emission (Fig. [31 top panel) shows a fairly 
regular geometry. About 15" north and south of the nucleus the molecular gas forms two gas 
lanes that run roughly in east- west direction over a length of 16" (north) and 20" (south). The 
southern lane is about twice as bright than the northern lane. In addition to these extended 
structures, several fainter and more compact clumps are detected. All major components are 
indicated in Fig. HJ The most prominent clump is closest to the nucleus (labeled as central 
clump in Fig. UJ) at a distance of 9" (~600pc) to the North. Four additional smaller clumps 
are present, three near the northern lane and the last one at about 25" west (western clump 
on FigSD of the center of NGC4151. 

We measured the CO fluxes within different components marked in Fig. H] using the tasks 
BLOT (to define the regions) and FLUX in GIPSY. All the fluxes presented here are corrected 
for primary beam attenuation. We then converted the CO(l-O) flux into molecular gas mass 
for each component, u sing a standard Galactic CO-to-H 2 conversion factor of Xco = 2 x 



10 20 cm 2 [K km s l ] 1 (iSolomon fc Barrettlll99ll ) and applying a mass correction for helium 



of 36% . The derived CO fluxes and molecular gas masses are listed in Tab. [3j In our entire 
map, the total flux recovered is 14.7 Jykms -1 and the corresponding cold molecular gas mass 
is about 4.3 xlO 7 M , with 60% of this mass being contained in the southern lane. Our PdB I 



data recover more flux than the JCMT single dish measurements (jRigopoulou et al.l 119971 ). 
They measure a total CO(2-l) flux of lco = 1.56 Kkms -1 , corresponding to Sco(i-o) 
I I .1 v km s 1 in integrated flux, assuming optically thick molecular gas, in the inner 23" 
(~1.5kpc). The total H 2 mass derived from the single dish data is then M(H 2 )= 2.4 x 10 7 M & 
which corresponds to a total molecular mass of M mo i = 3.2 x 10 7 M , after applying the mass 
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correction for helium. Therefore, the JCMT data missed about 25% of flux coming from the 
outer edge of the CO gas lanes. 

Interestingly, we detect no CO emission in the central 3", where warm H 2 emission i s 
observed via its NIR emission lines ( iFernandez et al.l Il999t IStorchi-Bergmann et al.l 120091 ) . 
We discuss the discrepancy between the CO and NIR H2 observations in the central regions 
in Sect. PI 



3.2.2. CO kinematics 

The CO(l-O) velocity field and the velocity dispersion map are presented in the middle 
and bottom panels of Fig. [3j The velocity field in the extended southern and northern 
lanes shows an ordered gradient along the lanes with values similar to those derived from 



HI observations (lAsif et al.lll998l ; iMundell et al.lll999h and an average velocity dispersion of 



approximately 7 kms -1 . Fig. [5] shows the integrated spectra of the three main components: 
the northern and southern gas lanes as well as the central clump. The integrated spectrum of 
the northern lane is centered on 80 km s -1 and the central clump at 55 km s -1 . The spectrum 
of the southern lane shows the presence of two kinematic components originating from each 
half of the lane - one centered on kms -1 , corresponding to the eastern part of the lane where 
it curves to the North, and the other centered on —40 kms -1 arising from its most prominent 
western part - separated by a region in the centre of the lane that shows an increased velocity 
dispersion, close to the kink point. Fig. [6] shows a selection of spectra taken close to this 
high velocity dispersion region. Asymmetry in the line profiles is clear, suggesting multiple 
kinematic components, but this region is spatially only marginally resolved while a second 
line component can be easily identified in several spectral. This boundary region between the 
two halves of the gas lane therefore may represent a region in which gas clouds with distinct 
velocities overlap along the line of sight or a genuinely mixed, disturbed zone. Larger velocity 
dispersion in the east part of the northern lane may also indicate that the eastern part of 
this lane is a distinct cloud and although the line profiles in this region deviate from a simple 
Gaussian shape they show no clear evidence of two components in this region. 

We extracted pv diagrams along three different positions that are presented in Fig. [71 
The pv diagram along the major kinematic axis of PA=22° shows strong emission from the 
gas lanes at +80 kms -1 and —50 kms -1 . The CO (1-0) line emission shows a velocity gradient 
indicating that the gas is participating in circular motion. Interestin gly this grad i ent is not 



symmetric with respect to the systemic velocity of 995 kms" 1 from iPedlar et al.l (119921 ) . it 
presents an offset of +15 kms -1 . In order to investigate the gas motion within the gas lanes 
two slits were placed along them (number 2 and 3, Fig. [7]). While the velocity gradient is very 
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low in the northern lane with a mean velocity of about +80 kms -1 , two distinct features are 
present in the southern lane. The eastern part of the lane shows constant velocities consistent 
with the systemic velocity, the western part exhibits a clear velocity gradient from -20 to 
-60kms _1 . 

The CO(l-O) rotation curve of NGC4151 was derived from the velocity field using the 
task ROTCUR in GIPSY. This task fits the velocity, assumed to be circular, within tilted 
rings. We run ROTCUR for rings of 2" widths from 3" to 19" radial distances centered on 
the AGN as traced by the peak of the 3mm continuum. The region defined by an angle 
of 5° around the minor axis was excluded for the fit. The inc lination and position angle 



were fixed at i=21° and PA=22° as derived from the Hi data bv Mundell et all (1199a ). The 
systemic velocity, v sys , was allowed to vary in a first ste p and found t o be at 1010 km s" 1 , 



i.e. offset by about +15kms _1 from the value derived by iPedlar et al.l (ll992f ). This offset is 
consistent with the offset in the velocity gradient along the major axis, as discussed above. 
Indeed, with a systemic velocity of lOlOkms -1 , the pv diagram along the major kinematic 
axis (PA=22°) presents a more symmetric velocity gradient, the southern and northern CO 
lanes lying then at -65 and + 65kms _1 , respectively. 

The resulting rotation curve, extracted with v sys = 1010 km s -1 is presented in Fig. |HJ 
The point at 9" presents a velocity lower than expected for circular rotation. This point 
corresponds to the north-eastern part of the southern lane, which has been shown to have 
kinematics distinct from the west part of the lane (Fig. The low velo city of the first point 



at 3" is due to beam smearing. Comparison to the Hi rotation curve ( IMundell et al.lll999[ ) 



shows a good agreement between the CO and Hi rotation curves (Fig. [BJ right). 

Therefore the overall CO kinematics are generally consistent with rotation in the plane 
of the galaxy. To demonstrate this further, Fig. M shows velocity as a function of position 
angle measured directly along the CO arcs. The sinusoidal curves overplotted show the 
expected variation of velocity along a circular annulus in the plane of the galaxy (i=21°, 
PA=22°) at a radius of 18" corresponding to the semi-major axis of the dust/CO ellipse. 
Although not exact, the correspondence between the plotted curves is good, ruling out the 
interpretation that the ellipticity traced by the CO arcs is a result of gas lying in a highly 
inclined circular disk tilted out of the main galactic plane as well as strong non-circular 
motions within the disk. The morphology is therefore intrinsically elliptical in the plane of 
the galaxy, as s uggested by the analysis of the Hi velocities along th e dust arcs (Fig.3 o f 



Asif et al.lll998l ) and consistent with predictions of gas flows in bars (jAthanassoulal Il992l ) . 
The rotation curve in Fig. [8] is thus only an approximation to the true kinematics and the 
offset between the derived CO and HI systemic velocities is a likely indication of disturbed 
CO gaseous kinematics in these arcs. 



-8 - 



Finally, the velocities within the central clump present larger deviations from circular 
motions, up to 50kms _1 below the circular velocity at a radius of 5"(Fig. [91 right). The 
central clump shows redshifted kinematics (see Fig. [3D significantly inconsistent with those 
expected for circular motion in this region, close to the minor kinematic axis of the galaxy. 
These observed deviations are suggestive of streaming motions, however, the uncertainties 
are too large to derive a more quantitative measure. Very deep CO imaging in the future 
is required to confirm whether this clump is an isolated cloud or part of a more continuous 
fluid flow towards the nucleus. 



4. Discussion 



4.1. Role of the Large-Scale Stellar Bar 



In order to analyse the importance of the large-scale stellar bar or oval for fueling the 
AGN, we need to establish the link betwe en the bar and the molec ular gas observed in the 
central 2.5 kpc. The detailed Hi study of iMundell fc Shone! (119991 ) showed that the stellar 
bar/oval is unusually gas-rich and that its Hi gas dynamics are consistent with the presence 
of X2 orbits in an oval bar: sharp velocity changes across the bright regions close to the 
leading edge s of the bar are fou nd which can be explained by offset shocks predicted by 
simulations (lAthanassoulal Il992l ) and provide dire ct evidence for the presence of these X\ 
and X2 families of stellar orbits. The Ha; image (IKnapen et al.l 12004 ) reveals Hll regions 
along the leading sides of the large-scale oval coincident with the peaks in the HI surface 
density. These H II regio ns are associated with the ionized and neutral gas in the streaming 
shocks ( lAsif et al.l 120051) that c o nnect to the CO gas lanes to form an inner gaseous spiral 
(Fig. [TUJ left). IMundell et al.l (119991 ) derived the resonance curves from the Hi velocity 
field and estimated the Inner Lindblad Resonance (ILR) to be at 2.8±0.6kpc. The ILR 
(red circle in Fig. [TUl left) coincides with the Hll regions, while the CO lanes are well 
inside the resonance. This picture is consistent with the molecular gas tracing a spiral 
density wave, driven by the la rge scale bar and extending from the ILR down to the center 
( lEnglmaier fc Shlosmanl l2000h . As we do not detect CO all the way down to the nucleus, 
this may argue in favor of th e presence of a second ILR th at limits the gaseous wave to exist 
between the two resonances ( lEnglmaier &: Shlosmanll2000f ). However, the sparse distribution 
of the CO gas, the lack of H I in the central regions and the non-linear behaviour of gas in 
a barred potential prevent useful contraints on the presence and location of an inner ILR 
from being derived directly from the gaseous rotation curves. Full dynamical modelling and 
numerical simulations constrained by millimetre, radio and optical velocity fields and K-band 
derived mass models will be presented in a future paper (Dumas et al. in prep.) and may 
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provide insight into location of an inner ILR and its effect on the gas dynamics. 



The CO gas lanes also coincide with two dust lanes seen in a V-I color map (jAsif et al. 



19981 ). as shown in Fig. [TDl right. The two red arcs deline ate a 11" x 18" dust ellipse 
inside the large oval/bar and also coincide with Hi emission (jAsif et al.l 119981 ) . The peak 
column densities in the Hi arcs are (0.7 — 1.8) x 10 21 cm -2 implying an atomic gas mass of 
Mhi ~ (0.13 — 3.3) x 10 6 M . Thus the ratio of atomic to molecular gas is 0.3 in the northern 
lane and 0.006 in the southern lane, which is in agreement with the val ues of the atomic-to- 
molecular gas rat io observed in the central regions of spiral galaxies (ISchuster et al.l 120071 ; 
Leroy et al.ll2008l ). with the southern lane being at the lower end. Moreover, the kinematics 
of the CO gas lanes and the H I arcs are in good agreement as shown in the previous section. 
Therefore the morphology and kinematics of the CO northern and southern lanes suggest 
that the molcular gas lanes are the continuation of the atomic and dusty ring, along elliptical 
orbits of the stellar bar. 



It is interesting to note that while the structure of the gas lanes is very well matched, no 
obvious counter-structure is present in the V-I map for the nuclear clump. High resolution 
HST/ACS imaging shows the presenc e of dust in the central 1.5 kpc, but no clear structure 
is seen in the dust (IBentz et al.ll2006l ). Connecting the gas flows in the large scale bar with 
the nuclear dust structures therefore remains a challenge. 



4.2. Cold and warm molecular gas in the innermost region 



In the central 2'.'5, excited H 2 line emission centered on the nuc leus and extended along 
the minor axis of the galaxy (Q','8 ~ 50 pc) ha s been observed by iFernandez et al.l (119991 ) 
and more recently by IStorchi-Bergmann et al.l (120091 ). The later showed that H 2 is ther- 
mally excited either by X-rays from the AGN or by shocks along the accretion flow to- 
wards the nucleus. This central H 2 emission has been interpreted as tracing the reservoir of 
(cold) gas from which th e supe r-massive black hole in NGC4151 is being fed. In particular, 
Storchi-Bergmann et al.l ( 120091 ) found a mass of hot H 2 gas of 240 Mm and based o n the 
conversion of warm-to-cold H2 mass ranging between 10~ 7 and 10~ 5 ( iDale et al.ll2005l ). they 
estimated the total mass of molecular gas (hot plus cold) to be in the range of 2.4 xlO 7 - 
1O 9 M . We should easily have detected such an amount of molecular gas. Indeed, the la 
upper limit for the CO(1-0) emission within the central 3" is 0.033 Jykms -1 , which translates 
into a limit for the cold molecular gas mass of 10 5 M , well below the lower limit expected 
from the NIR H 2 line. 



This discrepancy between warm and cold molecular gas mass can be interpreted in two 
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opposite ways. On one hand, it indicates that no significant amount of cold molecular gas is 
present close to the nucleus. This implies that the warm molecular gas observed there is not a 
good tracer for the general amount of (cold) molecular gas present. Thus the use of a general- 
ized ratio of warm-to-cold H 2 could be misleading in t he vicinity of AGN, which woul d explain 
the large over-estimate of cold molecular gas mass by lStorchi-Bergmann et al.l (120091 ) . On the 
other hand, in X-rays dominated regions (XDR), equilibrium molecular abundances may be 
affected by the X-rays radiation from the AGN. This leads in particular to a depleti on of CO 
abundance with respect to other molecules such as HCN (IBoger fc Sternberg 120051 ). Indeed 
many Seyfert galaxies show enhanced c e ntral HCN(l-O) emission such as, e.g . NGC1068 



(Jackson et al. 



1993; Tacconi et al 



1994; 



He 



(IKohno et all 120031 ). NGC5194 (IKohno et al 



fer fc Blitz 1995; Usero et al. 2004), NGC1097 



1999 ). NGC6951 (IKrips et alJl2007ah . In par- 



ticular, the latter study showed that the CO line emission largely underestimates the amount 
of cold molecular gas in the central region of the Seyfert 2 NGC6951. Therefore, CO itself 
may be a poor tracer of cold molecular gas close to an active nucleus. Under this assumption, 
the change of chemistry in the central kpc because of the AGN itself could explain the lack 
of CO emission inside the dust/gas lanes and the excitation of the NIR H2 would be then 
driven by the AGN rather than the accretion flow. Further study of molecular abundances 
in the central 3" of NGC4151 are needed to quantify the cold molecular mass in the central 
region of NGC4151 and to explain this discrepancy between (cold) molecular gas masses 
inferred from NIR H2 and CO line emission. 



4.3. Implication for AGN fueling 



The presence of non-circular motions in the central clump indicates that fueling of the 
AGN or at least a replenishment of the gas reservoir closer to the AGN might be ongoing. 



Anal ysis of near-IR images reveal no evidence for a second inner stellar bar ( IQnken et al. 



20071 ). This suggests that the disturbances of the regular motions might be caused by a 
transient eve nt, such as an i nterac tion with a nearby galaxy. Such an interaction has been 
proposed by iMundell et al.l (119991 ) to explain the distortion of the large two-arm stellar 
spiral. NGC4156 and UGC07188 have small projected distances to NGC4151 and there- 
fore may be candidate companions for such a galaxy interaction, although at respectively 
=6750 km s -1 and 6850 km s -1 the systemic velocities of NGC4156 and UGC07188 are 



J sys~ 



significantly larger than that of NGC4151. 

Our findings are in agreement with recent results from the NUGA survey that found no 
direct i ndication for ongoing fueling within the central hundred of parsecs for most of their 



targets (jGarcia-Burillo et al.ll2005t iBoone et al.ll2007l ; lLindt-Krieg et al.ll2008l ; ICasasola et al. 
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20081 ; lHaan et al.ll2009l ). It is, however, interesting that most NUGA galaxie s have a molecu- 
lar gas reservoir within the inner kiloparsec with mass b etween 10 7 - 10 10 M ^ (iGarciVBurillo et al 



2003b: Combes et al 



2004 



Garcia-Burillo et al 



2005 



Kr ips et al. 



2005 



Boone et al 



2007 



Krips et al.ll2007bt ICasasola et al 



2008 



Garcia-Burillo et al 



2009 



Combes et al 



20091 ). It 



is also interesting to note that most NUGA targets are either Seyfert type 2 or LINER 
galaxies. While gravitational mechanisms, driven e.g. by large scale stellar bars or galaxy 
interaction, are efficie nt in transporting gas to t he central kpc of galaxies, other processes 
(e.g. viscous torques, iGarcia-Burillo et al.ll2005l ; lHaan et al.l 120091 ) must take over to bring 



the gas closer to the nucleus and fuel the AGN. Moreover, recent nuclear star formation 



and AGN activity have been shown to be linked (IGonzalez Delgado et al.l 1200 ll ; iRiffel et al. 
20071 ; iDavies et al.l 120071 ). and stell ar outflows, from starburst on tens of parcsec scales, can 



also support the nuclear activity (IDavies et al.l 120071 ; iHicks et al.l 12009; Schar tmann et al. 



oi pj 



2009 ). NGC4151 presents no evidence for a recent nuclear starburst (jsturm et al. 19991 ). In 



fact, the star formation rate acr oss the whole ga laxy is unusually low given the presence of 
significant quantities of Hi gas (lAsif et al.ll2005j). therefore this sce nario appears to be un- 
likely for this particular Seyfert galaxy. Finally, iDumas et al.l ( 120071 ) showed that the ionised 
gaseous kinematics of Seyferts are more disturbed with respect to the stellar kinematics in 
the central kiloparsec than those of inactive galaxies, i.e. that kinematic disturbance scales 
with accretion rate in the inner kiloparsec where dynamical and activity timescales become 
comparable. 



5. Summary and Conclusions 

We have presented the first interferometric map of the molecular gas reservoir in NGC 4151 
using PdBI. A total molecular gas mass of 4.3 x 1O 7 M is detected in the central 2.5 kpc. 
Most of the gas is located in two gas lanes at lkpc distance from the nucleus. A nuclear 
clump containing about 3% of the total mass is found 600 pc north of the nucleus. However, 
no CO emission associated with the AGN has been detected down to a limit of 10 5 M for the 
central 3" indicating that estimates of the cold gas reservoir based on NIR H 2 line emission, 
based on the assumption of the ro- vibrational H 2 emission line in thermal equilibrium, can 
lead to severe overestimation of the amount of total gas present. 

The morphology and kinematics of the gas lanes are consistent with being driven by 
the large-scale stellar bar and also coincide with dust lanes seen in optical color maps. The 
detection of non-circular inward motion in the nuclear clump suggests that this cloud may 
trace gas that will ultimately flow towards the SMBH. However, the kinematics of the central 
clump could not be linked directly to those induced by the large-scale stellar bar. Detailed 
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dynamical modeling of the gravitational potential of the bar and the inner kpc regions is 
required to combine the molecular and ionised gas and stellar components in a coherent 
framework and identify the mechanism responsible for the observed non-circular motions, 
which is beyond the scope of this paper. Thus we speculate here that the potential inflow is 
caused by a transient phenomenon and the fueling process of the AGN might be erratic and 
intermittent. In agreement with recent studies of other nearby AGN, these findings suggest 
that no single (gravitational) mechanism is responsible for the fueling of nuclear activity. 
Finally, in the central kpc where dynamical and activity timescales become comparable, a 
fine balance may be required between fueling and feed back, but in NGC4151, which has 
an estimated accretion rate of m = 1.3 x 10~ 2 M Q yr -1 ( IStorchi-Bergmann et al.l 120101 ). the 
upper limit on cold molecular gas in the central 3" (10 5 M ) given by our CO observations 
is still sufficient to support the current level of activity of NGC4151 over an AGN lifetime 
of 10 7 yr- 
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Fig. 1. — Channel maps of the naturally weighted CO(l-O) line emission in NGC4151. The 
channels are 6.5 km s _1 wide and contours are plotted at -6, -3, 3, 6, 12 and 24<r with la 
= 2mJybeam _1 . The velocity marking in the top left corner of each panel is relative to the 
observed central velocity of v sys = 995 kms -1 . The CLEAN beam of 3.44" x 2.97" is shown 
in the bottom left corner of each panel. The cross marks the phase center of the observation. 
The dirty beam and the uv coverage are presented in the last two panels. 
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Fig. 2. — 3mm continuum map of NGC4151 using uniform weighting. The beam size of 
2.8"x2.2" is indicated in the bottom-left corner. The contour levels run from 1 to 12 mJy 
beam -1 with a step of 1 mJy beam -1 (3.5cr to 43cr in steps of 3.5cr).The cross marks the 
phase center of the observations, while the peak of the continuum coincides with the location 
of the AGN, within our astrometric uncertainty. 
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Fig. 3. — Moment maps of the CO(l-O) line emission in NGC4151: integrated intensity 
(top), velocity field (middle) and velocity dispersion (bottom). The contour levels are 0.25, 
0.5, 0.9, 1.3, 1.5 and 2 Jy beam" 1 kms" 1 in the intensity map; from -50 kms" 1 to 80 kms" 1 
in steps of 5 kms" 1 for the velocity field and 5, 7.5 10 and 14 kms" 1 in the dispersion map. 
The beam size (3.4" x 3.0") is shown in the bottom left corner of eac h panel and the cros s 
marks the location of the AGN as given by its 21cm radio continuum (IMundell et al.l 120031 ). 
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1 "" Northern lane 




Right Ascension (J3000) 



Fig. 4. — The individual components of the CO(l-O) distribution discussed in the text 
are indicate in the intensity map. The contour marks the 0.21 Jybeam -1 kms -1 level and 
roughly outlines the areas used to measure the CO(l-O) line flux for each component (see 
Table E} • 



- 21 - 




Fig. 5. — Integrated CO spectra of the three major components of the CO distribution: 
northern gas lane (left), southern gas lane (middle), and central clump (right). For all 
spectra, the x-axis is relative to the systemic velocity of v sys = 995 km s -1 (in kms -1 ) and 
the y-axis is the integrated flux density in Jybeam -1 . Note the second component in the 
profile of the southern gas lane that is centered at Okms^ 1 and mostly arising from the 
northeastern tip of this lane. 
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Fig. 6. — CO emission line profiles from different regions of the southern lane. The central 
panel shows the velocity dispersion accross the southern lane, between 3 and 15kms _1 from 
blue to white colors. The spectra have been extracted from single pixels of size 0.5" x 0.5", 
the locations of these pixels are indicated by the arrows. 



-23 - 




Fig. 7. — CO pv diagrams extracted along several position angles and slits in NGC4151. 
Contour levels go from —5a, — 4a, — 3a, 3a to 75a in steps of la = 2.8 mJy beam -1 . The 
locations of the slits are indicated in the intensity map (left) and the slits are labeled. The 
pv diagrams shown are along (1): the major kinematic axis (PA=22°); (2): the northern 
(PA=0°, offset=-7"; 14") and (3): southern (PA=55°, offset=-1.5"; -13") gas lanes (from left 
to right). 
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Fig. 8. — Rotation curve of NG C4151 derived from the CO(l-O) (black circles) and Hi 
velocity fields (gray crosses, from iMundell et al.l Il999l ). Both panels show the combined 
CO rotation curve (derived with v 



sys 



lOlOkms -1 ) and HI rotation curve, the left panel 
showing a zoom on the central region. The point at R=9" in CO rotation curve (left panel) 
corresponds to the north-east part of the southern lane, which has been shown to lie at the 
systemic velocity (Fig. [5]) and appears to be at a lower velocity than expected from circular 
rotation. 
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Fig. 9. — Left: Velocity along the CO lanes, as a function of position angle 6 around the lanes, 
from north to east. Red circles (respectively black squares) correspond to the northern lane 
(respectively southern lane). The two curved lines correspond to sinusoids of pure circular 
velocities for an inclined disk with PA=22° and i=21° at the major-axis (18") of the gas 
lanes. Right: Velocity within the central clump, as a function of position angle 9, from 
north to east (red points). The plain and dashed curves corresponds to sinusoids of pure 
circular velocities with disk orientation of PA=22° and i=21°, at a radius of 5", position 
of the central clump. In both panels, the sinusoids correspond respectively to a systemic 
velocity of v sys = 1010 kms -1 (solid blue line), and v sys = 995 kms -1 (dashed black line). 
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Fig. 10. — CO (1 -0) contours ove rlaid on the Ha map (left; iKnapen et al.ll2004j ) and a V-I 
color map (right: lAsif et al.lll998l ). The ellipse outlines the large-scale stellar bar present in 



NGC 4151 and the circle corresponds to the position of the ILR as derived by lMundell et al. 



(119991 ) . Enhanced Ha emission is present along the leading sides of the bar (indicated by 
arrows) that connect to the CO gas lanes. The V-I map was smoothed to 1.5" resolution to 
enhance the contrast for the dust lanes. Redder colors are presented by darker shading. The 
ENLR is evident as the very blue elongated shape in the nucleus. In both panel, the cross 
marks the position of the AGN. 
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Table 1. General Properties of NCG4151 



Property 


Value 


Reference 


Right ascension (J2000) 


12 fe 10 m 32. s 58 


Mundell et aL (2003) 


Declination (J2000) 


+39°24'21'.'06 


Mundell et al. (2003) 


Classification 


SAB(rs)ab 


RC3 fde Vaucouleurs et al. 1991) 


AGN Type 


Seyfert 1.5 


NED 


Distance 


13.3 Mpc 


Mundell et al. (2003) 


Scale 


1"=65 pc 




Systemic Velocity ~V S ys 


995 km s" 1 


Ulrich (2000) 


Inclination (i) 


21° 


Mundell et al. ( 1999) 


Position Angle (PA) 


22° 


Mundell et al. ( 1999) 


Bar Extent 


3'x2' 


Mundell et al. ( 1999) 


Bar PA 


130° 


Mundell et al. ( 1999) 
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Table 2. Dataset Parameters 



Datasct 


Weighting 


Channel Width 


Clean Beam 


RMS 






[km s~ 1 ] 


["x"] 


[mjy beam -1 ] 


continuum 


natural 


1150.5 


3.44x2.96 


0.20 


continuum 


uniform 


1150.5 


2.83x2.16 


0.28 


line-only 


natural 


6.5 


3.44x2.96 


2.0 


line-only 


uniform 


6.5 


2.83x2.16 


2.8 



Note. — Beam sizes and noise levels of the PdBI continuum maps and line-only 
data cubes derived at different velocity resolution and weightings. 
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Table 3. CO Fluxes and Molecular Gas Masses 



Component 


ScodV 
[Jy km s _1 ] 


M mo ; 
[1O 7 M ] 


Total 


14.7 


4.3 


Southern lane 


8.7 


2.5 


Northern lane 


5.0 


1.5 


Central clump 


0.4 


0.12 


Western Clump 


0.63 


0.18 



Note. — Integrated CO(1-0) line fluxes and 
corresponding molecular gas masses for differ- 
ent components identified in the CO distribution. 
The location and extent of the individual com- 
ponents are indicated in Fig. [4] The fluxes have 
been corrected for primary beam attenuation and 
the masses include 36% of Helium fraction. 



